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0^ . A simple numerical method to study the effect of an apphed magnetic field on the 

energy spectrum of non-periodic superlattice structures is presented. The magnetic 
O \ field could be either parallel or perpendicular to the growth direction. Our method 

\m) ' is based on the transfer matrix technique and on the effective mass approximation. 

CN ■ We discuss the advantages and disadvantages of the proposed approach using several 

. examples. In particular, we study the perturbation to the energy spectrum of periodic 

superlattice induced by the introduction of an enlarged well. We found that these 

J^ \ perturbations are negligible for B//z but relevant for B_Lz. Preliminary results for 

fS| ' Fibonacci superlattices in magnetic fields are presented as well. In these quasi- 

^? ■ periodic structures the energy levels become strongly dispersive in presence of a 

perpendicular magnetic field. 
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I. INTRODUCTION 



C^ 



X 

^ [ In 1970, Tsu and Esaki proposed superlattices as new structures with peculiar proper- 

ties (minibands) which allowed new physical phenomena to be studied. [0 In recent years, 
the improvement of epitaxial growth techniques renders disposable samples to study the 
properties of transport in the vertical direction (i. e. parallel to the growth direction, 
z). Resonant quantum tunneling through double barrier heterostructures, transport along 
superlattice minibands, and optical studies in biased multiquantum well structures have 
recently attracted great attention. [Q In short period superlattices, the coupling between 
adjacent wells leads to transport through superlattice miniband states. In addition, the 
present availability of high magnetic fields stimulated many experimental studies of mag- 
netotransport. [Q A problem of increasing interest is the effect of a transverse magnetic 
field, B, on the tunneling through a barrier separating either two semiconductors [Q or two 
superlattices 0. Magnetotransport in transverse fields has also been studied in a semiclas- 
sical macroscopic framework 0. Other topics concern the effect of an applied magnetic 
field perpendicular to the semiconductor layers of an heterostructure |§|J^. In this case, the 
magnetic quantization in the layer plane collapses the two-dimensional (2D) electron gas in 



a discrete set of ciclotronic orbits of zero degrees of freedom (OD system). 

The problem of the influence of an external magnetic field on the superlattice energy levels 
have been throughout studied by theoretical techniques. Two different approaches have been 



used: the first is based on a numerical integration either of the Schrodinger equation |T^ 



or of a reduced Luttinger Hamiltonian |Tl| and the second solves the Schrodinger equation 
using the method of expansion of the solutions in series of analytical functions (like sine 
functions [0, ipergeometrical confluent functions [^ or parabolic cylinder functions [0). 
In this paper we propose a different approach based on the extension of the well-known 
Transfer Matrix (TM) method ]T5| , p!B| to superlattices in magnetic field. The interest in our 
approach is the possibility to study superlattices formed by whatever a succession of layers. 
The TM method has been widely used to solve the problem of non-periodic structures, fl^ 



Its simplicity is related to the use of the product of only 2x2 matrices. This renders the 
method very appealing because it is easily implementable on small computers. 

In the Appendix we present the details of our approach to the Transfer Matrix method in 
absence of magnetic field. For simplicity we have restricted our discussion to the electronic 
problem, i. e. we neglect the holes. In section II we study the case of B\\z, where z 
is the growth direction: the superlattice miniband width is varied due to the magnetic 
quantization in the heterostructure layer planes, i.e., to the conservation of the Landau level 
index along the different layers. In section III, we analyzed the case of BJ^z. A new equation 
which couples the motion in the layer plane and in the vertical direction should be solved 
to obtain the energy levels. Then, we study two particular problems: the enlarged well 
in a superlattice (section IV) due to its importance for vertical transport studies and the 
Fibonacci superlattice (section VI) as an example of non-periodic superlattice. In section 
V we have computed the density of states for magnetic fields oriented either parallel or 
perpendicular to z. Finally, section VII ends the paper with a summary of the main results. 

II. MAGNETIC FIELD PARALLEL TO Z 

Through the effective-mass approximation it is possible to separate the electronic wave- 
function in a component (usually assumed to be a plane wave) for the motion in the layer 
plane and in a component (solution of the equation (|TT])) for the motion in the z direction. A 
magnetic field applied in the z direction does not introduce a mixing of the parallel and the 
perpendicular electron wave functions. However, it changes the plane waves for the motion 
in the layer plane into harmonic oscillator wave functions and it collapses the free electron 
dispersion relation in a discrete set of Landau levels. [^ In this case, the eigenvalues £"„ are 

En = E, + nuj^{n + 1/2) (1) 

with n = 0,1, 2, .. the Landau level number, Ez the eigenvalue for the motion in the vertical 
direction, u^ = eB/m*c the cyclotron frequency for the GaAs effective- mass m*, e the 
electronic charge and c the light velocity in vacuum. In (|I]), we have neglected the spin 
magnetic energy. 



To obtain E^, it is necessary to solve the effective- mass Schrodinger equation 



2m' 



■V' + Vr,{B,z) 



ij{z) = E^z) (2) 



where the magnetic field is included by changing the ordinary superlattice potential, Vsl{z), 
with an effective one-dimensional potential, Vn{B,z), different for each Landau level n, 

Vn{B, z) = Vsl{z) + h{n + 1/2) [u;^ - u;,] = K + Hn + 1/2) [cu^ - Ub] (3) 

in the Al^.Gai_^.As barriers and 

Vn{B,z) = VsL{z) = (4) 

in the GaAs wells. Vq is the conduction band discontinuity between GaAs and Al^Gai_^As, 
and Uh is the cyclotron frequency for the Alj.Gai_a;As effective-mass. The potential described 
in equations (^ and (^) reflects the conservation of the Landau level number through the 
layers, f^ Consequently, E^ depends on the applied magnetic field. 

Figure 1 shows the highest and lowest electron energy level of the miniband of a 30/30 
superlattice (30A thick wells and barriers, with x = 0.33 in the barriers) versus the applied 
magnetic field for different values of n. [|l^ By increasing the applied magnetic field, the 
electron energy levels move to low energy. The energy shift is larger for higher Landau 
level number. This effect is explained by equation (^). The Al^Gai_xAs effective mass is 
higher than the GaAs effective mass, so that the cyclotron frequency in the Al2;Gai_a:As 
layers is lower than in the GaAs layers. Thus, the second part on the right hand side of 
equation (|^) is a negative term proportional to B and n. li B or n increases, the potential 
barrier is lowered, and the electron energy levels are shifted to lower energies. Hence, 
the coupling of the electron wave functions in adjacent wells is augmented increasing the 
miniband width. To evaluate this effect, we compare the dependence of Ez on B for various 
n with the corresponding ciclotronic kinetic energy hujy^{n + 1/2). For B measured in Tesla, 
huw = 1.72 B meV. So that the conservation of the quantum number n changes Enhj 5 %. 
This effect could be increased if a higher aluminum concentration is used. 

III. MAGNETIC FIELD PERPENDICULAR TO Z 

In this case, the Hamiltonian describes a spinless particle in a layer of the superlat- 
tice, with an effective mass m* and charge e, subject to a constant and uniform transverse 
magnetic field B, i.e. 

where p = —ikV, and A is the magnetic vector potential. We fix the direction B\\x to 
simplify the equations. The Hamiltonian takes the form 

/i^ 2 e^B^z^ _ iehB^d_ 

2m* 2171*0^ m*c dy 

where the gauge A = (0, —Bz, 0) has been used. By using a plane wave in the x and y 
direction and neglecting spin effects, it is possible to find the envelope wave- function tp{z) 
that describes the motion of the carriers along the z direction of the heterostructure through 
the equation 



H^{z) = E^{z) (7) 

with E the electron energy and the effective mass Hamiltonian, H, given by 

where the definitions ^o = —hkyc/eB have been used, ky is the wave-vector in the y-direction. 
Note that the eigenvalues E depend on zq. 

To calculate the electron energy eigenvalues for different zo (dispersion relation) we have 
assumed a flat band approximation for the potential Veff{z). This consists in taking Veff{z) 
constant in each layer. For a layer i between two interfaces Zj and %+i, we chose a mean 
value for the magnetic potential. We define Zi = {zj + Zj+i)/2 and the second term on the 
right hand side of equation (H) yields a constant value Vi = m*u'^{zi — zo)'^/2. Then, in 
order to use the transfer matrix method, we substitute E^E — V^ in each formula of the 
Appendix. This approximation is based on the large number of barriers and wells that forms 
a superlattice. The averaged fiat band potential is expected to be a good approximation 
of the real potential Veff{z) if the wave function extends over many wells and barriers. In 
this case, we can neglect the parabolic shape of the potential Veff{z) in each individual 
well and barrier. Increasing the number of layers, the accuracy of the averaged fiat band 
approximation is increased. In this way, the application of a parabolic magnetic field on a 
superlattice is modelized by means of a variable height fiat band potential. This method is 
similar to the approximations that R. Tsu and L. Esaki have used to compute the effect of 
an external electric field on a superlattice. ||20| 



Our approximation has a low infiuence for low magnetic fields. In addition, we have found 
a very small numerical differences for the dispersion relations of electrons in superlattice 
under high magnetic fields (up to B=20T) with respect to the results of J.C Maan []TU[. 



In Ref. [jlO[ the Schrodinger equation (P) has been numerically integrated. Moreover, to 
check further the approximation used in our calculations we have tried to reproduce more 
accurately the potential Veff{z) in each layers. In a barrier or well placed between the Zj 
and Zj^i interfaces we discretize the layer in a finite number (L) of regions. The interfaces 
between the new L regions are found at ze = Zj + (-Zj+i — Zj)i/L. We define a fiat band 
potential in each i region as Vf = Vsl{,Z() + m*uj'^{zi — zqY /2. Increasing the number L, 
the new step- like potential will be closer and closer to Veff{z). Table I reports the electron 
energy levels in a 30/30 Alo.asGao.erAs-GaAs superlattice with B = 5T for different L values 
and a fixed zq. Including a large number of quantum wells in the superlattice (20 wells) the 
energy levels obtained with L = 2,3, ... are slightly different from the results with L = 1. 
The L = 1 case represents a very good compromise between computational efficiency and 
the accuracy of the results. For L ^ 1 the number of the matrix products are changed from 
A^ — 1 (N is the number of layers) to (A^ — 1)L. 

We have plotted in Fig. 2 and 3 the calculated electron dispersion relations {zq versus E) 
using 20 wells and L = 1 for a 30/30 and 60/60 superlattices at S = 5T and 4T, respectively. 
These are representative of two electronic dispersion relations with one miniband (30/30 
superlattice) or two minibands (60/60 superlattice). The dispersion relation is flat if the 
eigenvalues belong to the minibands calculated in absence of magnetic fleld and it is almost 
parabolic if not. |T^ Anti-crossing in the energy levels occurs for several Zo values. For the 
60/60 superlattice, it is found that the anticrossing occurs also between the energy levels 



of different minibands. From a semiclassical point of view, the anticrossings between both 
minibands correspond to cyclotron orbits of electron states localized in different minibands. 
In this case, one electron could tunnel from one to the other superlattice minibands. 

In the literature, other methods based on the transfer matrix method have been proposed 
to study the effect of the transverse magnetic field on the superlattice energy spectrum (see 
e. g. |ir3| , p^ ). However, these other methods use some analytical solutions of equation (|]) in 
each layer (ipergeometrical confluent functions [|T^] or parabolic cylinder functions |T^). The 



numerical calculation of these functions requires long CPU times and introduces numerical 
errors and instabilities that in our scheme are avoided. Using the analytical method of Ref. 
||r3| to calculate the energy levels for a double barrier structure one needs 10 times longer 
CPU time than to compute the same structure with the present method with L = 6. 

IV. THE ENLARGED WELL PROBLEM IN SUPERLATTICES 



Recently, Chomette et al |^ have proposed to use an enlarged well (EW) as a marker 
of the arrival of carriers at a certain distance from the surfaces. This method allows to 
study vertical transport in superlattices by means of optical techniques. Including a well 
in the superlattice, which is wider than the superlattice wells, introduces a localized state 
in the energy spectrum. The level due to the EW has a lower energy than the superlattice 
miniband. In this way, the photoluminescence from the EW can be used as a measure 
of the transport through the superlattice. [^ In absence of any applied magnetic field, it 
is found that the original superlattice miniband is slightly affected by the introduction of 
an EW, so that, the results obtained with this technique can be used to understand the 
transport through superlattice minibands without an EW. The enlarged well technique has 



been used to study magnetotransport in superlattice in the cases of -B| |z |]22| and B^z. 
In this section, we analyze the perturbation of the superlattice energy spectrum due to the 
EW when an external magnetic field parallel and perpendicular to the growth direction is 
applied. 

For parallel magnetic field, the electron energy levels are easily found through the method 
of calculation described in section II. Table II reports, for B=50 T and n=0, the electronic 
energy levels for a 30/30 Alo.ssGao.erAs-GaAs superlattice and for the same superlattice with 
an EW of 60A in the middle. The electronic energy levels are only slightly changed by the 
presence of the EW even for this very high B value. Consequently, it is possible to use an 
EW as a marker to study magnetotransport in superlattices with B\\z by means of optical 
techniques. 

For B^-Z and using the method of section III, we show in Fig. 4 and 5 the energy levels 
for the same superlattices as those used to obtain Fig. 2 and 3 except that a 60A and 120A 
EW are introduced in the middle of the 30/30 and 60/60 superlattice, respectively. The 
introduction of the EW breaks the periodic symmetry of the superlattice and introduces a 
weaker inversion symmetry respect to the center of the EW. The interaction between the 
unperturbed dispersion relation of the superlattice and of the EW results in a situation very 
different from that shown in Fig 2 and 3. In particular, the flat dispersion reminiscent of 
the electronic miniband becomes dispersive after the introduction of the EW. This is due 
to the anticrossing between the EW states and the SL states. This perturbation of the SL 
energy levels disappears far away from the EW region. In a semiclassical point of view. 



the anticrossing between the EW energy level and the superlattice miniband corresponds 
to skipping orbits centered in the barriers on the left and right sides of the EW. [^ One 
electron in a superlattice state can tunnel to the EW localized state. As a consequence of this 
interaction, for B±z, the EW effect on the superlattice energy spectrum should be considered 
when the EW technique is used to study magnetotransport through a superlattice. 

V. DENSITY OF STATES 

In quantum wells, the two dimensional density of states, p{E), is usually taken as a step 



function. ||2^ In a superlattice composed by a finite number of layers, p{E) is a series of 
step, one for each level in the miniband. A parallel magnetic field collapses the 2D energy 
levels in a discrete set of OD Landau states. In this case, the step like function for p{E) is 
changed to a discrete set of Dirac (5-functions. p4| Including a Gaussian broadening of the 
5-functions caused from the presence of disorder, p{E) results 

p{E) = CY^j:^exp{-[E^ + huUn + 1/2) - E]ya'} (9) 

where C is a normalization factor, E^ are the energy levels of the superlattice in the z 
direction with the index m running over the superlattice miniband, n is the landau level 
index and a is the gaussian width. We have neglected the second term on the right hand 
side of equation (|) to obtain the energy eigenvalues Em, i.e. to compute equation (||) we 
have used the E^. values obtained in absence of a magnetic field. Figure 6 reports p{E) for 
a 30/30 superlattice for different magnetic fields. The Landau quantization changes p{E) 
in a set of isolated peaks. For each energy level in the miniband, the distance between two 
Landau states is given by the cyclotron frequency, huw Increasing the applied magnetic 
field, the cyclotron frequency, which is proportional to B, increases, and then, the distance 
between different peaks. 

For B±z, the density of states is. 



P^^^ = ^797^ ^ r \ h\F % f.^1 ^'-^ - E^{^o)]dzo (10) 

c (2vrft)^^Jo \i 2\E - Em[zQ)\ 



where Em,{zo) are the solutions of equation (|T^) which depend on the cyclotron orbit center 
Zq, and La;, L^ are the superlattice dimension. The square root dependence of p(E) is due 
to the ID motion in the x-direction which is affected neither by the superlattice potential 
nor by the magnetic field. A singularity in the denominator of the square root occurs at 
E = Em{zo) and it leads to important modifications in the optical and transport properties 
of these systems. 

Figure 7 shows different p(E) computed using eq. (|10D- For a 60/60 superlattice at B=4 
T (full line in Fig. 7), the existence of two electronic minibands is observed at about 60 
meV and 220 meV. The second electronic minibands is wider than the first one and this 
allows the resolution of distinct peaks in p(E). These peaks show a characteristic l/yE 
dependency. In Fig. 7, we compare also the p(E) results for B=5T of a 30/30 superlattice 
in presence (solid line) and in absence (dashed line) of a 60 A EW. In the miniband energy 
region (~150 meV), the effect of the EW induces a loosing of the ID typical divergence, 



especially for the states near the high energy edge. The direct contribution of two non- 
dispersive EW eigenvalues in the density of states is marked by two arrows. If the B value is 
increased no more non-dispersive levels exist and the ID divergence in the miniband region 
will disappears. 

VI. THE FIBONACCI SUPERLATTICES 

Another interesting problem which can be tackled by the present approach is quasi- 
periodic superlattices. These structures do not exhibit translational invariance but long 
range order. The interest stems from the fact that the Bloch theorem is inapplicable and 
that they provide an intermediate class between periodic superlattice with extended elec- 
tronic states and randomly disordered superlattice with exponentially localized electronic 
states. Merlin et al. |2^ proposed to grow quasi periodic superlattice using the Fibonacci 
sequence. These superlattices, called Fibonacci superlattices (FSL), are interesting because 
using a simple recurrent law one can obtain series of increasing period starting from usual 
simple superlattice and ending with a completely non periodic superlattice. [^^-^| We re- 



call that a quasi-periodic FSL of order / is constructed by / applications of the Fibonacci 
transformation law A -^ AB and B ^ A starting from A, yielding the following sequences: 
B (/ = -1), A (/ = 0), AB (/ = 1), ABA (/ = 2), ABAAB (/ = 3), ... . 

In Fig. 8 we show the electronic energy dispersion for different FSL and for B = and B 
= 5 T. The evolution with the magnetic field strength of the energy level scheme is shown in 
Fig. 9 for a FSL with index f=l and for B = 0, 1, 3, 5 T. In our calculations, following Ref. 
^,TB], the block A was formed by 28 A of GaAs and 30.8 A monolayers of Alo.23Gao.77As, 



and the block B by 28 A monolayers of GaAs and 67.2 monolayers of Alo.23Gao.77As. We 
have repeated the Fibonacci series to obtain a total superlattice thickness of ~ lOOOA for 
the different FSL but the f=5 FSL which has a superlattice thickness of 1500 A. 



For B=0, our method yields the same results as the calculations of Ref. ||27|jl6|l . In 
summary, two kinds of states are formed in the spectrum of a FSL. The first is formed by 
single-A pattern (typical of / = 1 FSL), the second kind is formed by double-A pattern 
(typical of / = 2 FSL). A double-A pattern is analogous of two coupled QW separated by a 
thick barrier. Little overlap exists between these two kinds of states for the / = 1 and the 
/ = 2 FSL, whereas the electronic energy spectrum of the other FSL presents overlaps with 
either the / = 1 FSL energy bands or / = 2 FSL energy bands (see in Fig. 8 the B=0 T 
results for f=3 or f=5). The effect of an external magnetic field is the breaking of this simple 
scheme. In particular, for /=0 (periodic superlattice) and B= 5 T we get non- dispersive 
states in the region of the superlattice miniband and dispersive states for higher energies. 
While for / 7^ 0, almost all the states are dispersive. The quasi-periodicity of the FSL is 
lost and the energy dispersions are very similar for different /-values. The magnetic field 
localizes the electronic orbit in finite spatial regions. Hence the electron are subject to the 
local potential (K//)- K// does not show the quasi periodicity of the Fibonacci superlattice 
potential. This fact is shown in Fig. 9 for /=!. By decreasing the electron cyclotronic 
orbit (i. e. increasing B) increases the dispersion of the states. For B=l T the energy 
dispersion is flat. Increasing the magnetic field the energy dispersion changes: for B=3 T 
we observe both dispersive and non dispersive states while for B=5 T almost all the states 
are dispersive. We note that the band of states at k, 250-300 meV which develops for B= 5 



T and for f=0-3 in Fig. 8 is due to the finite size of the superlattices. 

VII. SUMMARY 

In summary, we have presented a simple method to compute the effect of an external 
magnetic field on the energy spectrum of superlattices. This method is based on the transfer 
matrix approach. Due to the strong non-parabolicity and the valence band mixing of hole 
states in superlattices, we have applied our approach only to the electronic states in this 
paper. 

A magnetic field parallel to z leads to a full quantization of the electronic motion in the 
layer plane. The conservation of the Landau level indexes along the superlattice increases 
the energy miniband width and induces a dependence of the energy levels on the Landau 
index and the applied magnetic field value. The density of states is changed from a step-like 
function to a set of (5-like functions. 

For B^z case, we have computed the dispersion relations and the density of states for 
electrons in a superlattice. The energy spectrum shows the existence of energy regions with 
dispersive and non dispersive energy states. The density of states is characterized by a 
typical 1/yE dependency for non-dispersive states and by a smooth increase for dispersive 
states. 

We have applied our method to study the enlarged well problem in superlattices with 
an applied magnetic field. With B\\z, we have found that the presence of the EW is only a 
small perturbation on the superlattice energy spectrum which permits the use an EW in the 
study of magnetotransport through superlattices. For B^.z case, the EW strongly perturbs 
the superlattice energy spectrum. 

Fibonacci superlattices which are a particular class of non-periodic structures show some 
interesting aspects due to the loss of the quasi-periodicity in magnetic fields. Our preliminary 
results suggest that due to the increase of the magnetic localization of the electrons the self- 
similar energy spectrum typical of these quasi-periodic structure is lost. A magnetic field 
of 5 T is enough to render dispersive all the electronic states while for periodic superlattice 
the energy dispersion is still flat at this value of the magnetic fied. 
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TRANSFER MATRIX METHOD 

Through the effective-mass approach it is possible to separate the wave function in a 
component parallel to the layers and a component in the growth direction, z. The eigenvalues 
in the z direction, Ez-, and the envelope wave functions, ip{z), in a superlattice can be 
calculated using spherical masses and parabolic dispersions in the layer plane. PU| , P(]| In this 
way, the motion of electrons and holes in the growth direction can be decoupled and reduced 
to the solution of a one- dimensional Schrodinger equation for each kind of particle. Let us 
assume to solve the problem of electrons. The Schrodinger equation is: 

8 



h' 



2m* (z) 



V^ + Vsl{z) 



tpiz) = Eip{z) 



(11) 



where m* is the effective mass for electrons in the different semiconductor layers and Vsl{z) 
is the superlattice potential. The superlattice potential is given by the difference of the 
energy gap between different layers which we have approximated using a one dimensional 
Kronig-Penney-like potential. [^ The energy zero is defined at the bottom of the GaAs 
conduction band for the electrons. The envelope wave function solution of (|11]) in the zth 
well is 



ilji{z) = AiCos{kiz) 



mj 
Bi— — sin[kiz) 



(12) 



where Ai and Bi are two constants, ki = {2m*EY''^ /h and m* is the effective mass in the 
ith layer. The envelope wave function in a zth barrier is 



*'^' - Tr" ^ w^^'- 



(13) 



where k, = J2Tn*{Vo — E)/h, with Vq as the height of the potential barrier. 

Using a carrier current conserving boundary condition, the continuity of ip{z) and 
ifj'{z)/m*, the coefficients Ai and Bi of the layer i are joined to those of the layer z+1 
by a transfer 2x2 matrix Si^i^i of determinant 1 defined in the following way |1^ 




S,. 



n+l 




(14) 



Si^i+i is equal to {Cf\^i)^^Ci if the layer z is a barrier and the (z-/-i)th layer a well, and to 
{Cij^i)^^Cf if the ith layer is a well and the {i+l)th layer a barrier. Both functions C are 
calculated at the interface between the i and (i+l) layer. C^ is defined as 



C 



cos{kiz) ^sin{kiz) 



Ki 



—^sin{kiz) cos{kiz) 



(15) 



and Ci as 



Cf 



(16) 



By consecutive applications of the matrix iSj^j+i we form the matrix Ai which joins the 
first and last layers of the superlattice 




-(«;)^n/-(«: 



(17) 



with N the number of superlattice layers. Imposing vanishing wavefunctions in the first and 
last superlattice layers, we find that Ai = B^ = 0, and from (p!7|) the (2,2) element of the 
matrix Ai results in 



9 



M22 = (18) 

The solutions of the eigenvalue equation ([T8|) yield the miniband spectrum for electrons in 
the superlattice. 

This method, which use the product of 2x2 matrices, can be easily implemented on 
small computers and it permits to calculate any finite succession of layers and not only 
periodic superlattices. On the other hand, we neglect energy band non-parabolicity for the 
electrons and the heavy-hole light-hole mixing in the valence band. Consequently a big error 
is expected when these effects are important. To estimate the errors on our calculations, we 
have compared our results with the ones of Ref. 1^^ where energy band non-parabolicity 
and valence band mixing have been considered. The differences between the results for short 
period superlattices (period shorter than 120 A) are of about 2% for electron energy levels 
and of about 10% for hole energy levels. Due to the neglection of valence band mixing, we 
have found a worst agreement for the heavy hole energies than for electron energies, where 
the agreement is rather satisfying. To go beyond our approximations, larger matrices are 



necessary, |32] which result in a larger computational effort. 



10 
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FIGURES 

FIG. 1. Lowest and highest electronic energy level in the miniband of a 30/30 (30 A wide 
wells and barriers) Alo.ssGao.eyAs-GaAs superlattice versus the applied parallel magnetic field B 
for different Landau level index. Solid line: n=0, Dotted line: n=l, Dashed line: n=2. The zero 
of the energy scale is the bottom of the GaAs conduction band. 

FIG. 2. Solid line: electron dispersion relation for the 30/30 (30 A wide wells and barriers) 
Alo.ssGao.eyAs-GaAs superlattice for a perpendicular magnetic field B=5T. Dashed line: position 
of the equivalent superlattice potential. The zero of the energy scale is the bottom of the GaAs 
conduction band. 

FIG. 3. Solid line: electron dispersion relation for the 60/60 (60 A wide wells and barriers) 
Alo.ssGao.eyAs-GaAs superlattice for a perpendicular magnetic field B=4T. Dashed line: position 
of the equivalent superlattice potential. The zero of the energy scale is the bottom of the GaAs 
conduction band. 

FIG. 4. Solid line: electron dispersion relation for the 30/30 (30 A wide wells and barriers) 
Alo.asGao.eyAs-GaAs superlattice with a 60A enlarged well in the middle, for a perpendicular 
magnetic field B=4T. Dashed line: position of the equivalent superlattice potential. The zero of 
the energy scale is the bottom of the GaAs conduction band. 

FIG. 5. Solid line: electron dispersion relation for the 60/60 (60 A wide wells and barriers) 
Alo.ssGao.eyAs-GaAs superlattice with a 120A enlarged well in the middle, for a perpendicular 
magnetic field B=6T. Dashed line: position of the equivalent superlattice potential. The zero of 
the energy scale is the bottom of the GaAs conduction band. 

FIG. 6. Density of states in arbitrary units versus energy for a 30/30 (30 A wide wells and 
barriers) Alo.asGao.erAs-GaAs superlattice in a parallel magnetic field. Solid line: B=1T. Dashed 
line: B=2.5T.The zero of the energy scale is the bottom of the GaAs conduction band. 

FIG. 7. Density of states versus energy for a 60/60 (60 A wide wells and barriers) 
Alo.ssGao.eyAs-GaAs superlattice (solid line) at B=4T. The density of states for a 30/30 super- 
lattice is also shown for B=5T without (solid line) and with (dashed line) an enlarged well in the 
middle. The magnetic field is perpendicular to the growth axis. The zero of the energy scale is the 
bottom of the GaAs conduction band. The arrows indicate the enlarged well bound states. 

FIG. 8. Electronic energy levels of different Fibonacci superlattices with index / as a function 
of the cyclotronic orbit center (zq) for a perpendicular magnetic field B= T (discs) and B=5 T 
(points). For a definition of the Fibonacci index / and for the superlattice parameters see the text. 
For B= T the energy levels do not depend on Zq. The energy barrier in the superlattice is 177 
meV and the zero of the energy scale is the bottom of the GaAs conduction band. 
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FIG. 9. Electronic energy levels of a Fibonacci superlattice with index /=1 as a function of 
the cyclotronic orbit center Zq and for different values of the perpendicular magnetic field B. For 
B= T the energy levels do not depend on Zq- For a definition of the Fibonacci index / and for 
the superlattice parameters see the text. The energy barrier in the superlattice is 177 meV and is 
indicated in the figure by an orizontal line. The zero of the energy scale is the bottom of the GaAs 
conduction band. 
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TABLES 

TABLE L First seven electron eigenvalues for 20 coupled quantum wells (30/30 superlattice) 
with B=5T at different L values. The cyclotron orbit center value has been taken in the middle of 
the sample. The zero of the energy scale is the bottom of the GaAs conduction band. 

Electron energy levels (meV) 

L^l 106.84 114.11 121.26 128.29 135.15 141.81 148.18 

L=2 104.13 113.20 119.32 126.27 134.13 139.24 145.28 

L=3 103.93 113.02 118.92 125.69 133.55 138.45 144.93 



TABLE IL First seven superlattice eigenvalues for 20 coupled quantum wells (30/30 superlat- 
tice) with B\\z and for the same superlattice potential with an implemented EW. The magnetic 
field is B=50T and Landau level index is n=O.The zero of the energy scale is the bottom of the 
GaAs conduction band. 



Electron energy levels (meV) 



SL 97.97 98.67 99.83 101.46 103.52 106.03 108.95 

EW 96.23 97.32 98.51 99.62 102.12 105.81 107.41 
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